Anomalous Quantum Interference Induced by Landau-Zener Transitions in a Strongly 

Driven rf-SQUID Qubit 
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We irradiated an rf-SQUID qubit with large-amplitude and high frequency electromagnetic field. 
Population transitions between macroscopic distinctive quantum states due to Landau-Zener transi- 
tions at energy-level avoided crossings were observed. The qubit population on the excited states as 
a function of flux detuning and microwave power exhibits interference patterns. Some novel features 
are found in the interference and a model based on rate equations can well address the features. 
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Coherent quantum dynamics of superconducting 
Josephson qubits driven by weak external held have been 
extensively investigated in experiments*!* 2 -'^'^ These 
experiments used Rabi oscillation to probe the macro- 
scopic quantum coherence of such systems in the time do- 
main, where the driving frequency equals to the energy- 
level separation. On the other hand, quantum coher- 
ence can also be investigated with a large-amplitude 
field, which can drive the qubit throughout its energy- 
level spectrum. In this case, Landau-Zener (LZ) transi- 
tion process occurs around each level-crossing position 
with a finite probability while the quantum evolution 
is adiabatic away from the level-crossing positions. 7,8 
Repetition of LZ transitions gives rise to Stuckelberg or 
Ramsey-type oscillation o 9 ' 10 in analogy to Mach-Zehndcr 
(MZ) interferometer. The MZ-type interference is a 
unique quantum coherence signature complementary to 
Rabi oscillation. 

Recently, coherent dynamics of superconducting qubits 
in the strongly driven regime dominated by LZ transi- 
tions were extensively studied, especially in two level sys- 
tems, providing an alternative method to manipulate and 
characterize the qubiti 11 ' 12 ! 13 ' 14 ! 15 A recent work makes 
a step forward^ : the superconducting flux qubit as a 
multi-level system was driven through several avoided 
crossings between energy levels. The qubit popula- 
tion under large-amplitude fields exhibited a series of 
diamond-like interference patterns in the space parame- 
terized by static flux detuning and microwave amplitude. 
Spectroscopic information can be obtained from the field 
amplitude dependence of the qubit population^ That 
experiment was done with a superconducting persistent- 
current qubit, which has relatively a small loop size 
and a long coherence time. The applied microwave fre- 
quency was much smaller than the energy-level spacing 
between discrete energy levels and therefore the resonant 
peaks/dips tend to merge into a continuous band. 

Comparing with a superconducting persistent-current 
qubit, an rf-SQUID qubit has a much larger geometric 
loop size (about 100 times larger in loop area) which 
usually leads to a very short coherence time (~1 ns). 



In order to observe LZ interference in such a large sys- 
tem, the time interval between subsequent LZ tunneling 
events should be less than the relevant coherence time, in- 
dicating a high-frequency field is required. In a previous 
work^ we observed the quantum interference fringes in 
an rf-SQUID qubit due to LZ transitions at one energy- 
level avoided crossing. In this paper, we report the ob- 
servation of simultaneous presence of two-set of interfer- 
ence fringes associating with LZ transitions at two nearby 
avoided crossings. In addition, the qubit population can 
be modulated by the field amplitude at certain static flux 
bias points. We developed a model based on rate equa- 
tions and addressed these features very well. 

Our sample was fabricated with Nb/AlO^/Nb trilayer 
on an oxidized Si wafer by using the standard pho- 
tolithography method. As shown in Fig. 1(a) and 1(b), 
the rf-SQUID qubit is essentially a superconducting loop 
with a second order gradiometric configuration inter- 
rupted by a small dc-SQUID. When an external magnetic 
flux <E>* colse to $o/2, where $o is the flux quantum, 
is applied to the qubit loop, the system potential takes 
a double well shape. The flux states in different wells, 
serving as the qubit states, correspond to macroscopic 
circulating currents flowing in opposite directions. Two 
on-chip current bias lines are used to control the mag- 
netic flux applied to the qubit loop (3>p and the small 
dc-SQUID loop ($^ JJ ) respectively. determines the 
asymmetry of the potential and $9 J modulates the po- 
tential barrier height. By adjusting and $y? JJ , the 
shape of the system potential can be fully controlled in 
situ. An additional on-chip current bias line is used to 
supply microwave (MW) pulses. The qubit flux states 
can be read out by a hysteretic dc-SQUID inductively 
coupled to the qubit loop. We can easily bias the readout 
dc-SQUID at its maximum sensitivity region by applying 
an external magnetic flux (<3?j c ). 

The sample was placed on a chip carrier enclosed in 
a superconducting aluminum sample cell and the experi- 
ments were performed in a dilution refrigerator at a base 
temperature about 20 mK. The device was magnetically 
shielded and all electrical lines were carefully filtered and 
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FIG. 1: (color online) (a) Optical micrograph of an rf-SQUID 
qubit (big loop interrupted by a small dc-SQUID with two 
Josephson junctions), (b) Schematic of the on-chip circuits 
(crosses represent the Josephson junctions), (c) The mea- 
surement time sequence . 



equipped with proper cold attenuators to minimize cir- 
cuit noise (see Ref^ for detailed description). 

In our experiment, we biased at a certain value 

such that the critical current of the small dc-SQUID is 
close to its maximum. Therefore the barrier between 
two potential well maintained at a relatively high value. 
There are several flux states below the top of the barrier 
in each potential well. The inter-well flux state transition 
rate between the ground states of two potential wells is 
very small at flux biases close to $o/2- Since the intra- 
well relaxation time is very short, the state of the qubit 
freezes quickly after turning off microwave irradiation. 
The time sequence of measurement is shown in Fig. 1 (c) . 
For each measurement cycle, a 10 fis microwave pulse is 
applied to the qubit loop to induce sinusoidal excursions 
through the energy levels about a static flux detuning 
(5<I>* = - $o/2. After a delay time about 1 fxs, the read- 
out is performed by driving the dc-SQUID with a current 
pulse comprising a 15 ns sampling current followed by a 
15 [is holding current just above the retrapping current of 
the dc-SQUID. By properly adjusting the magnitude of 
the readout current pulse, the dc-SQUID either switches 
to finite-voltage state or stays at zero voltage state cor- 
responding to the qubit states being in left well \L) or 
in right well \R). The measurement repetition frequency 
is 5 KHz. By repeating the trials for 4xl0 4 times, we 
obtained the average switching probability, representing 
the population in \L) state. By changing the flux detun- 
ing and microwave power step by step, we then measured 
the dependence of qubit population in \L) state on flux 
detuning and microwave power. 



Fig. 2(a) shows the resonant peaks and dips when ap- 
plying a microwave with frequency 17.0 GHz and nominal 
power -16 dBm. The peaks (dips) reflect the mcrowave 
induced excitation from flux state \R) to \L) (\L) to \R)). 
The width of the resonant peaks/dips is generally on the 
order of several m<E>o and the decoherence time of our 
system is estimated to be less than 1 ns. The strong 
population inversion is evident, suggesting that higher 
energy levels other than the first two lowest states are 
involved^ Without microwave the qubit step center is 
shifted away from the symmetry position of the energy 
spectrum (marked by the left dotted line in Fig. 2(a)) 
due to the presence of the circulating current in the read- 
out dc-SQUID^ Then we fixed the flux bias and swept 
microwave power. Fig. 2(b) shows the qubit population 
in \L) state as a function of the nominal microwave power 
at a static flux detuning = 3 m$o (marked by the 
right dotted line in Fig. 2(a)). The population exhibits 
a damped oscillation around 0.5. For two level systems, 
a good Bessel function dependence of the population on 
microwave amplitude was observed in a superconduct- 
ing persistent current qubit and the modulation period 
is proportional to the ratio of microwave amplitude and 
the microwave frequency, 1 — Therefore, in our experiments 
a large microwave power is required to observe the peri- 
odic behavior since the microwave frequency is very large. 
However, for a very high microwave power, the LZ transi- 
tions are not the dominant processes anymore. To make 
a more quantitative understanding, we hereby propose 
a simple model, which is illustrated in Fig. 2(c). The 
qubit is initialized in \0R), the ground state of right well. 
When the microwave amplitude is not very large (i.e., the 
double well potential always exits), the LZ transitions 
(with rate W) between \0R) and \nL) (the nth excited 
state of left well) , and the fast intra- well relaxation (with 
rate 7) from \nL) to |0L) dominate the dynamics of the 
qubit population. These processes cause qubit popula- 
tion inversion and lead to the non-monotonic dependence 
of qubit population on microwave power. When the mi- 
crowave amplitude goes up, it is much easier for the the 
qubits to be directly excited from the ground state to a 
higher energy level \N) (with rate g), which is above the 
top of the barrier and has a wavefunction spanning both 
potential wells. The qubit population in \N) tend to fast 
relax back to the ground states (\0L) and \0R)) of both 
wells with an approximately equal relaxation rate T and 
therefore this process equalizes the qubit population in 
the two wells. Taking all transition processes into consid- 
eration, one can obtain the rate equations to describe the 
qubit level occupations pi(i = 0, 1, 2, 3, corresponding to 
\0R), \0L), \nL) and \N) respectively) as: 

Po = -(W + g + T Q1 )p Q + Tiopi + Wp 2 + Tp 3 , 
pi = T Q1 p Q - (g + T 10 )p 1 + 7p 2 + I> 3 , 
P2 = Wp - (W + 7)pz, 

P0+Pl+P2+P3 = l (1) 

Under the assumption of classical noise and using per- 
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Nominal MW Power (dBm) MW Amplitude A (ui) 

FIG. 2: (color online) (a) Normalized qubit population in \L) state as a function of static flux detuning Four peaks and 
dips are induced by a fO fis microwave pulse at 17.0 GHz and -16dBm. As a guide to eyes, the left dotted line marked the 
midpoint of the peak/dip positions. Without microwave, the qubit transition step center is displaced from this midpoint. The 
right dotted line corresponds to static flux detuning 5$^ = 3m<I>o (b) Normalized qubit population in \L) state versus the 
nominal microwave power at a fixed flux detuning 3 m«l>o. (c) Schematic of the relevant transition processes involved in Fig. 
2 (b). \0R) and \nL) are the ground states of right well and left well respectively. \N) represents a higher energy level which 
is not localized in either potential well. The red (blue) lines show the excitation (relaxation) path. The purple lines represent 
the LZ transitions between \0R) and \nL). (d) Calculated qubit population in \L) state as a function of microwave amplitude. 
The relevant parameters used in our theoretical model are: u>/2n = 17 GHz, m = 8, A/2-7T = 7 MHz, F2/27T = 2 GHz, *y/2n = 
r/27r = 2 GHz, ri /27T = 0.6 KHz, r i/27r = 0.1 KHz, a/27r = 5 Hz, b = 1.4. 



turbation theory, the LZ transition rate between \0R) 
and \nL) takes the formil: 



W(e,x) 



r 2 Jl(x) 



2 ^ (e - tow) 2 + T 2 . 



(2) 



where A is the avoided crossing between \QR) and \nL), 
e is the dc energy detuning from the avoided crossing A 
and T2 — I/T2 is the system dephasing rate. J m (x) are 
Bessel functions of the first kind with the argument x = 
A/u>, where A and u> are field amplitude and frequency 
respectively. For a m-photon transition process, e = muj 
when on resonance, the LZ transition rate approximately 
reduces to: 

A 2 Jm{x) 



W(e,x) 



(3) 



The transition rate from the ground state to \N) is sen- 
sitive to the barrier height during the microwave driving 



process. We simply assumed the transition rate have an 
exponential dependence on field amplitude A, i.e., g = 
ae hA , where a, b are two fitting parameters. L10 (L01) is 
the slow inter-well relaxation rate between |0L) and \0R). 
In the stationary case, po = pi = P2 = i>3 = 0, Eq.(l) 
can be solved numerically with appropriate model param- 
eters. Fig. 2(d) shows the calculated qubit population of 
state \L) as a function of the microwave amplitude. The 
theoretical curve clearly exhibits the main features of the 
experimental results, indicating our model captured the 
underlying physics of the system. 

Fig. 3(a) shows the contour plot of the qubit pop- 
ulation in state \L) versus the positive static flux de- 
tuning (5<I>* and the nominal microwave power with a 
fixed microwave frequency at 17.0 GHz. According to 
the previous worksi^, we expect that the resonant peaks 
are equally spaced on the flux axis and the distance be- 
tween the adjacent peaks corresponds to the microwave 
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FIG. 3: (color online) (a) Dependence of qubit population in 
\L) state on static flux detuning <5<E>y and nominal microwave 
power, with microwave frequency at 17.0 GHz. Part of two 
overlapped interference sets can be observed and population 
inversion is clearly demonstrated, (b) Energy-level diagram 
illustrating the multi-photon transition processes. The qubit 
is initially in \0R), the ground state in right well. Microwave 
harmonically drives the qubit through the energy levels that 
are approximately linear in flux detuning and LZ transitions 
occur at the avoided level crossings Aor,uL and AoR,(n+l)L 
(another two level crossings locating at positive flux detuning 
are not shown in this figure). The red and green path, with 
different energy- level slope, represent the main processes to 
generate the first and second interference set respectively, (b) 
Numerical calculation of qubit population in \L) state versus 
flux detuning and microwave amplitude. The calculated inter- 
ference patterns show good agreement with the experimental 
observations (Fig. 3(a)). 



different energy-spectrum slopes, indicating two different 
transition paths are present. The two interference sets 
have overlap along microwave power axis. In addition, 
the distance of the resonant peak increases with the flux 
detuning. Therefore, the interference pattern looks a lit- 
tle complicated and is different from the regular patterns 
induced by Landau-Zener transitio n 11 ! 12 ! 13 ! 14 ' 15 ! 18 . In or- 
der to explain the interference we have to deal with the 
multi-level system. Since the barrier of the rf-SQUID 
is very high, the Landau-Zener transition rate is very 
low for the lower energy states. We thereby consider the 
avoided crossings close to the top of the barrier, which 
is qualitatively illustrated in Fig. 3(b). The qubit is 
initially in state \0R) at a certain static flux bias 
The first set is associated with LZ transitions at level 
crossing Ao_R iri L (also Aql^h, not shown in Fig. 3(a) for 
simplicity) and the red path represents the main transi- 
tion processes for generating the first set. The second 
set is associated with LZ transitions at level crossing 
A fl,(n+i)L (also A oij ( n+1)fl ) and the green path rep- 
resents the main transition processes for generating the 
second set. The fours peaks in the first set correspond to 
ni-photon and up to (ni+3)-photon transitions and the 
three peaks in the second set correspond to n2-photon 
and up to (n 2 +2)-photon transitions. The slope of the 
green path is smaller than that of the red path. There- 
fore the flux spacing of the resonant peaks of the second 
set is larger than that of the first set. On the other hand, 
the slopes of the energy levels close to the top of the 
barrier are flattened, leading to the increase of the dis- 
tance of the resonant peaks. Moreover, the most notable 
feature is the overlap of these two sets of the resonant 
peaks, which can also be understood from Fig. 3(b). For 
a high-frequency microwave, the peak region of the in- 
terference fringe of the first set can cover a wide range 
of microwave power. Therefore, during the modulation 
of the first set, the second level crossing can be easily 
reached and then the second interference set starts to 
show up. This feature is different from the observation 
of well-resolved individual interference diamonds when 
using a relatively low microwave frequency. 



frequency. With the amplitude increasing, the resonant 
peaks emerge one by one, starting from the flux close to 
$o/2. The flux position for each peak almost does not 
change while the population follows the Bessel functions 
of the first kind. However, here we observed two sets 
of the resonant peaks. The first set includes four res- 
onant peaks, locating at around flux detuning 3 m<I>o, 
12 m$ , 23 m$ , an d 37 m$ , respectively. The sec- 
ond set lies in the higher microwave power region and 
includes three resonant peaks, locating at around flux 
detuning 10 m<I>o, 22 m<I>o, an d 38 m<I>o, respectively. 
We can group them into two sets because as microwave 
frequency varies, their flux bias positions move with two 



We can quantitatively compare the experimental re- 
sults with the theory by performing numerical calcula- 
tion based on the rate equation approach introduced in 
Refi 2 ^. In our calculation, the six energy levels shown in 
Fig. 3(b) are taken into account. For convenience, we use 
|0), |1), |2), |3), |4>, |5) to denote flux states \0R), \0L), 
\nR), \nL), \(n + l)R), \(n+l)L) respectively. A 0R . nL 
and A jj („ +1 ) i are also denoted by Ai and A 2 . Then we 
can write down the rate equations to describe the qubit 
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level occupations pi(i — 0, 1, 2, 3, 4, 5) as 

Pa = -(Wo3 + Was + r i)p + TioPi + T(p 2 + pi) 

+ W 30 p a + W 50 p 5 , 
Pi = ToiPo - (Win + Wu + T w )pi + T(p 3 + p 5 ) 

+ IU21P2 + W41P4, 
P2 = W12P1 - (W 2 1 + T)P2, 

P3 = W 3P0 - (W30 + r)p 3 , 
Pi = Wupi - (W41 + T)pi, 

P0+P1+P2+ P3 + Pi + P5 = 1 (4) 

where Wij are the LZ transitions between states |i) and 
\j), which can be explicitly written as: 

w 03 = w 30 = w 12 = w 21 = ^Y, i ^#r^ - 

2 ' fei — mo; r + J- 9 

m / z 

2 ^— ' 69 — mcJr + 1 9 
(5) 

where ei and €2 are the dc energy detuning from the 
avoided crossing Ai and A2 respectively. T is the fast 
intra-well relaxation rate from the excited state to the 
ground state and Tio (T01) is the slow inter- well relax- 
ation rate between |1) and |0). In the stationary case, 
Po = pi = pi = P3 = pi = P5 = 0. In our calculation, 
uj/2tt = 17 GHz, n = 4, n x = 8, n 2 = 10, A 0RAL /2n 
= Ai/2tt = 7 MHz, Aok,5l/27t = A 2 /27r = 13 MHz. 
As shown in Fig. 3(b), these model parameters describes 
the energy-level structure, which can be determined from 
the parameters: qubit loop inductance L « 1.3 nH, junc- 
tion capacitance C « 35 fF and junction critical current 
I c w 610 nA. Other relevant parameters in calculation 
are: r 2 /27r = 2 GHz, T/27T = 2 GHz, and r i0 = 0.5 
KHzxe^ 10 (eio is the energy difference between |1) and 
|0}). By using the above parameters and taking into ac- 
count the different energy-level slope for different transi- 



tion paths, one can numerically solve Eq.(4) in the sta- 
tionary case and obtains Fig. 3(c), showing the qubit 
population in \L) state as a function of the microwave 
amplitude and static flux detuning. The agreement be- 
tween experimental results and numerical calculation is 
remarkable. 

It has been known that the quantum interference pat- 
terns in a superconducting flux qubit is dominated by 
the microwave frequency and qubit decoherence time^S. 
Due to the larger size, the rf-SQUID qubit usually has 
much shorter decoherence time than that of the supercon- 
ducting persistent current qubit. Therefore, in order to 
resolve interference pattern one has to apply microwave 
with large frequency. Since the response of the microwave 
is proportional to A/ui, a large microwave power is sub- 
sequently required. The large power and frequency will 
lead excitation to many energy levels. Although the con- 
tribution from high levels complicates the interference 
pattern, we can still quantitatively understand the inter- 
ference using LZ transition. For more theoretical discus- 
sions on the interference features with large-amplitude 
and high-frequency field, one may refer to our recent 
note.— 

In conclusion, we observed the quantum interference 
fringes in a superconducting rf-SQUID qubit driven by 
large-amplitude and high-frequency microwave field. The 
interference pattern exhibits two sets of overlapped reso- 
nant peaks with well-resolved multi-photon transitions. 
The two interference sets correspond to two different 
transition paths and each path is associated with LZ 
transitions at a particular level-crossing position. By con- 
sidering the energy levels close and above the top of the 
energy barrier between the potential wells, we numer- 
ically calculated interference patterns. The agreement 
between the simulation and experimental results reflects 
the multi-level feature of the rf-SQUID qubit. 
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